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Abstract: Heavy quarkonium at finite temperature has been the subject of intense theo- 
retical studies, for it provides a potentially clean probe of the quark-gluon plasma. Recent 
studies have made use of effective field theories to exploit in a systematic manner the hi- 
erarchy of energy scales that characterize the system. In the case of a quarkonium in a 
medium whose temperature is smaller than the typical momentum transfer in the bound 
state but larger than its energy, the suitable effective field theory is pNRQCDjjTL; where 
degrees of freedom with energy or momentum larger than the binding energy have been 
integrated out. Thermal effects are expected to break Poincare invariance, which, at zero 
temperature, manifests itself in a set of exact relations between the matching coefficients of 
the effective field theory. In the paper, we evaluate the leading-order thermal corrections 
to the spin-orbit potentials of pNRQCDjjrpj^ and show that Poincare invariance is indeed 
violated. 
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1. Introduction 



Heavy quarkonium can be systematically studied by means of non-relativistic effective field 
theories (EFTs) Recently, the non-relativistic EFT framework has been extended to 
allow the study of heavy quarkonium in a thermal bath [|, |, |, |, |, 0, |, |. The relevance 
of heavy quarkonium, as a probe of the hot medium created by heavy ion collisions at 
modern accelerator machines, has been remarked since long time [p^ . 

Heavy quarkonium is characterized by the scales typical of a non-relativistic bound 
state and the thermal bath is characterized by the thermodynamical scales. The former 
are the heavy-quark mass, m, the inverse of the typical radius of the system, 1/r, and the 
binding energy, whereas the latter are the temperature T (or rather multiples of vrT) and 
possibly smaller scales. The non-relativistic scales are hierarchically ordered: the mass is 
much larger than the inverse of the typical radius of the system, which in turn is much 
larger than the binding energy. In the weak-coupling regime, which may be relevant for 



the lowest quarkonium resonances |11, 12 1, the inverse of the typical radius of the system 
scales like mos, while the binding energy scales like ma^. In the following, we will assume 
that the quarkonium is a weakly coupled bound state and that the temperature is such 
that 

m » mas ^ T » ma^ ■ (1-1) 

Moreover, we will assume that all other thermodynamical scales as well as the typical 
hadronic scale, Aqcd, are either of the same order as or smaller than the binding energy. 
As it has been argued in [0, 13 1, this situation may be relevant for the phenomenology of 



the ground states of bottomonium in heavy-ion collisions at the LHC. A detailed study of 
the spectrum and the width of quarkonium under a more restrictive hierarchy than ( |1 . 1| ) 
can be found in 

One may exploit the hierarchy (|l.l| ) by systematically integrating out degrees of free- 
dom associated with the scales m, mog and T and by substituting QCD with low-energy 
EFTs better suited to describe the quarkonium in the assumed regime. It turns out that. 



under the scale hierarchy (1.1), quarkonium behaves like a Coulombic bound state with 
small thermal corrections induced by the thermal bath. 

After integrating out the scales m and mag from QCD, one obtains potential Non- 



Relativistic QCD (pNRQCD) u4, 15 1. Although in a non-relativistic EFT Poincare invari- 



ance is not apparent, still such invar iance is realized non-linear ly through a set of relations 



among its matching coefficients |16, 17, 18, |T^, 20|. One of these relations is the so-called 
Gromes relation that relates the spin-orbit potential with the static potential. 

Integrating out the temperature T from pNRQCD, one obtains pNRQCDjjrp^ [^, p2| . 
This EFT has the same degrees of freedom as pNRQCD, but the matching coefficients 
and, in particular, the potentials get thermal corrections, while the Yang-Mills Lagrangian 



modifies into the Hard Thermal Loop (HTL) Lagrangian [23|. Because the thermal bath 
introduces a preferred reference frame (the one in which the thermal bath is at rest), 
one expects that pNRQCDy-pj^ breaks Poincare invariance;^ this would have consequences 



""^More precisely, the thermal bath breaks Lorentz invariance. Since, however, in the context of non- 
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for the properties of a heavy quarkonium moving with a certain velocity relative to the 
thermal bath. In order to make the statement more quantitative, we calculate in the paper 
the leading-order thermal contributions to the spin-orbit potential of pNRQCDjjrpL 
check them against the Gromes relation. 

The paper is organized as follows. In Sec. ^, we write the pNRQCD Lagrangian and 
review how Poincare invariance is realized in the effective field theory at T = 0. Section |3| is 
devoted to the computation of the leading thermal corrections to the singlet static potential 
and to the part of the spin-orbit potential of pNEQCD^jr^^ that depends on the centre- 
of-mass momentum. In Sec. we show that the thermal corrections violate the Gromes 
relation. In Sec. |5|, we calculate the leading thermal corrections to the complete spin-orbit 
potential and, in Sec. we present our conclusions. 



2. pNRQCD 



The largest scale in the hierarchy (1.1) is the heavy-quark mass. Integrating it out leads 
to NRQCD 1 25, 26 1. The matching is insensitive to the lower-energy scales, such as the 
thermal ones, therefore the EFT Lagrangian is the same as the one at zero temperature. 

The next relevant scale is the inverse of the typical radius of the system, which, in the 
weak-coupling regime, scales like mas- By integrating out this scale, we obtain weakly- 
coupled pNRQCD [0, whose La grangian reads 



^pNRQCD = j d'r Tr {st {id^ -K)S + Qt {iDo - K) O 
{S^hsoO + H.C.) + C.C.j - [o^/iooO + C.C. 

O^hiOh^, + C.C.] } - + q^ ^Jpq: ■ 



(2.1) 



i=l 



The fields S = S lc/\/iVc and O = T"/-v/Tp, are the quark-antiquark colour-singlet and 
colour-octet fields respectively, nj \s the number of light quark fields, A^'c = 3 is the 
number of colours, Tp = 1/2 and iDqO = id^O — qAqO + OgAQ. The trace is intended over 
colour and spin indices; C.C. stands for charge conjugation and H.C. stands for Hermitian 
conjugation. Quark-antiquark fields depend on the centre-of-mass coordinate, R, on the 
relative distance between the quark and the antiquark, r, and on time. Gluon fields depend 
only on R and on time; this is achieved by multipole expanding them in r. The operators 
hg and ho do not contain gluon fields, except for the case of ho in covariant derivatives, and 
may be interpreted as the colour-singlet and the colour-octet Hamiltonians respectively; 
terms contributing to hg and ho are ordered in powers of 1/m. The operators hgo, hoo, h^o 
and h^g contain gluon fields; terms contributing to them are ordered in powers of 1/m and 
r. We will detail these terms in the following. Since, again, the matching is insensitive to 



relativistic EFTs exact relations among the matching coefficients have been derived from the Poincare 
algebra ^ , an approach that may be traced back to , in the paper we will keep referring to the 
broader Poincare invariance. 

^We adopt here and in the following the notation of ]1S|]. 



- 2 - 



the lower-energy scales, the Lagrangian (|2.l| ) is the same as the one at zero temperature. 
We assume to be in the laboratory reference frame, which we define as the frame where an 
infinitely heavy quarkonium would be at rest. 
The singlet and octet Hamiltonians read 



(1) 



h. 



m 



4m 



+ 



y(2) 
Vs,o 



m 



+ ... , 



(2.2) 



where m is the heavy-quark mass, P = — iDp, and p = —i'Vr- Terms in the Hamiltonian 
have been ordered in powers of 1/m; the dots stand for higher-order terms. We recall 
that the non-relativistic power counting has 1/r ~ p 
centre-of-mass momentum, P, may be as large as T. 
The static potentials read 



mas and vj^J 



~ mQg, while the 



r IJSr r 



(2.3) 



where Cp = {N^ — 1)/ (2iVc) and ay^ and ay^ are series in a^, ay^ is known up to three loops 
1 27, 28 1, whereas ay^ up to two loops [^]. At leading order, it holds that ay^ = ay^ = Og. 
The singlet and octet propagators of pNRQCD can be expanded as 



^singlet 



E-hs-T.s{E)+ij] E-hs + iv 

i 



+ 



Qoctet 
'~>ab 



(E) 



E-ho- ^o{E) + ir] 



E -hs + i-q 



^s{E) 



1 



E - hs+ir] 



ab 



E -ho + irj 



+ 



E - ho + ir] 



^o{E) 



1 



ab 



E - ho + irj 



+ 



(2.4) 



(2.5) 



where S^^o are the colour-singlet and colour-octet self-energies. According to the power 
counting of pNRQCD and the fact that vj}J ~ vj^J <C mOg , the propagator 1/ [E—hs^o+if]) 
may, in turn, be expanded as 



E - hs,o + E- hfl + ii] 



E - hf}o + in 



.(2) 



4m 



+ 



m 



E - hf}o + iv 



(2.6) 



where h^}'^ = — Cp— and /i^"-* = ^ — I — — ^ are the leading-order singlet and octet 

m r m 2Nc r 

Hamiltonians respectively. 

The explicit form of the non-static potentials vj^^ and vj"^^ can be read from 

In particular, we will concern ourselves with the colour-singlet, Vlss, and the colour octet, 

(2) (2) 

VlSo, spin-orbit potentials, which are part of Vs and Vo respectively. They can be 
conveniently split into a part that depends on the centre-of-mass momentum P and a part 
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that depends on the relative momentum p: 

VlSs = ^ '—^^ '-VLSsa{r) + ^ ' \ ^ -VLSsh{r), (2.7) 

(r X P) • (crW -cr(2))^^ (r x p) • (crW + cr(2)) 

where the Pauh matrices cr^^^ and cr^^^ act on the heavy quark and antiquark respectively. 
At leading order, it holds that 

T/ I \ CpOs , - SCpOs , - 1 as 3 as 



2 r-^ 2 r'^ 4iVc r^* 4A'c r'^ 

(2.9) 

which implies that Vls sa ~ Vissb ~ V^soa ~ Visob ~ "i^af. 

The matching coefficients, and among them the potentials, that appear in the pN- 
RQCD Lagrangian obey a set of relations due to Poincare invariance For the kinetic 
terms in Eq. ( |2.2D , such relations impose that the coefficient of the operator P^/(4m) 
is one, which also implies that the coefficient of the operator p^/m is one. Among the 
relations fulfilled by the potentials, there are some exact relations linking the spin-orbit 
potentials Vls sa aiid Visoa to the static potentials: 

VLSsa{r) = , VLSoair) = "^^^ , (2.10) 

where /(r)' = df[r)/dr. The first relation is known as the Gromes relation, because it was 
first derived in ||2l| from the transformation properties of some Wilson loops under Lorentz 
boosts (see also ||l8[| ). 

In the next section, we will compute the leading thermal correction to the spin-orbit 
potential ViSsa by matching to pNRQCDhtl- To this end we will need some of the terms 
appearing in the operator hso of the pNRQCD Lagrangian ( |2.1| ). These may be ordered in 
powers of 1/m and r as 

hso = /.ir^ + + h^i^'^ + hg''^ + /.if) + . . . , (2.11) 

where the indices refer to the order in powers of 1/m and r respectively. The dots 
stand for higher-orders terms. The explicit expressions of h^so^ may be taken from |19] and 
read 

1/(0:1)/- 

hf^'^^ = -Y^^r . gE , (2.12) 

(r • -«) (r • gB) - LY,2±S^r • gB , (2.13) 
2m m 2r 

hg''^ = ^ Vjy\r) {p., r xgB} + ..., (2.14) 
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him 

1 1 an 



16m^ 



Fi„Y)(r)a«.[Px,5E] 



1 



+ 



16777,2 
1 



2 (r.*T«){P-,(5Exr)} 

(r) 



r 

T/(2,0) 



+ 



+ 



1 



^ sob"'\' ) 



1 



T/(2,0) 
v soe 



(r) 



8m2 



{(r-5E),P-(r x<T«)} 
-{(r •P),cr(^) • (r x c/E)} 
{P-,r X5B} + ... . 

,(0,2) 



(2.15) 



Charge conjugation invariance requires that hsd = 0. The field E is the chromoelectric 
field, E"^ = the field B is the chromomagnetic field, = —eijkF^''/2, [Px,(7E] = 
P X gE — gE x P and similarly for the anticommutators. For hi],'^^ and /lio'^^ only the 
P-dependent terms have been displayed. The coefficients cp and Cg are inherited from 
NRQCD and encode non-analytical contributions in 1/ m, whereas the various 
come from the matching to pNRQCD and encode non-analytical contributions in r. At 



leading order in the coupling, the matching gives cp = Cg 



Vsl''^\r) = 1, while all other matching coefficients are of order as or smaller. 
Poincare invariance imposes further constraints on the matching coefficients. For in- 
stance, at the NRQCD level, we have that 2cf — — 1 = to all orders |16, 17]. At the 



lB.ndvt^\r) = V^lf\r) 



level of pNRQCD, the following exact relations hold |19] 



,y(2,0)( 



Combining the last two it follows that 



(2.16) 
(2.17) 

(2.18) 
(2.19) 



An interesting consequence of this relation is that, since Vso'^^f) is at least of order a. 
1 31] but has not infrared divergences at that order [32|, V^JJ^??(?') is at least of order a^. 



3. pNRQCDhtl 

In this section, we compute the leading temperature-dependent correction to the colour- 
singlet static potential, V^^'^^ and spin-orbit potential, Vissa, which amounts to matching 
the colour-singlet static and spin-orbit potentials in the EFT that follows from pNRQCD 
after integrating out the temperature. This EFT is called pNRQCDjjrpj^ ^, 22]. Its 
Lagrangian reads 

/^pNRQCDhtl = / Tr {s^ [i^o - - 6Vs\ S + 0+ [iDo - K - dVo] o} 

+Chtl + ..., (3.1) 
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Figure 1: The leading heavy quarkonium self-energy diagram. Single lines stand for leading-order 
quark-antiquark colour-singlet propagators, ijiE ~ h^^'^ + ir]), double lines for leading-order quark- 
antiquark colour-octet propagators, iSab/{E — h^o^ + irf), curly lines for gluons and the vertices are 
the chromoelectric dipole vertices induced by h^so^\ 



where 6Vs and 5Vo are the thermal corrections to the colour-singlet and colour-octet po- 
tentials respectively, £htl is the HTL Lagrangian for the gauge and light-quark degrees of 



freedom [23] and the dots stand for singlet-octet and octet-octet operators suppressed by 
powers of 1/m and r, whose computation is beyond the scope of the paper. 

The terms 5Vs and 5Vo encode the corrections to the potentials induced by the ther- 
mal bath. They are evaluated by matching real-time propagators in pNRQCD with the 
corresponding expressions in pNEQCD^rp^. It has been pointed out in |^ that the 
doubling of degrees of freedom, typical of the real-time formalism, does not affect heavy 
quarks for which the unphysical degrees of freedom decouple. Hence, in the following, we 
will only deal with the physical singlet and octet propagators, which have been written in 
Eqs. (2^) and ( |2.5D respectively. The matching may be done perturbatively, because we 



have assumed T ^ Aqcd; contributions coming from other thermal scales do not affect 
the potentials, because we have assumed that the other thermal scales are at most as large 
as the binding energy. More specifically, the pNRQCD singlet propagator (|2.4|) is matched 
to the pNRQCDjjrpL expression 



E - hs - 5Vs + ir] E - hg + ir] 



There is no self-energy contribution in (|3.2| ), because this would correspond to a scaleless 

1 /2 

integral eventually irrelevant (e.g. in dimensional regularization it would vanish). Zs = 
1 -|- dZg is the normalization of the singlet field in pNRQCDjjTj^; SZg is of order Os and 
a function of r, which implies that it does not commute with hs- At our accuracy, 6Zs is 
real. 

We will assume the thermal bath to be a quark-gluon plasma at rest with respect to the 
laboratory reference frame. This implies, in particular, that the Bose~Einstein distribution, 

which describes the distribution of the gluons in the bath, depends only on their energy. 
3.1 Singlet static potential in pNRQCDjjrp^ 

We start by evaluating the leading thermal correction of the static potential, dVg'^'^ , which 
was first done in [S], and the leading thermal correction of 5Zs, which is new. These 
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corrections originate from the pNRQCD diagram shown in Fig. ||, whose contribution to 
T,s{E) reads 



2 2 



kQ — k'^ + if] 



d^k i 2 

(2vr)^ E - -ko + ^r^ ' 

+ 2tt5 {k^o-k^)nBi\ko\) 



(3.4) 



where k is the modulus of the three-momentum k. Evaluating the integral over the mo- 
mentum region ko, kr^T:$> (E-h'^o^) implies that we may expand 



E-h^^^ 



E - h^^^ 



ko + ir] 



+ i 



E-h'o 



(0) 



-ko + irj {-ko + irjy {-ko + irj)'^ 



+ 



(3.5) 



It is convenient then to regularize the integral in dimensional regularization, because the 
non-thermal part, which is scaleless, vanishes, and we are left with the thermal part only, 
which is finite. The leading contribution comes from the term in ( |3.5| ) that is linear in 
E — h^o^ ; after performing the integration, it reads 



A useful identity is 



27r 



(3.6) 



E-h^!>^ 



which follows from /i, 



1 

2 

(0) 



+ 



\{r^E-hf^}-{v(')-VP)r\ (3.7) 



(0) 



(0) 



v.. 



(0)n 



The identity is 



useful, because the first term in the right-hand side is of order 1/m and, therefore, does 
not contribute to the static potential, the second term contributes only to the singlet 
normalization, Zs, and the third term to the static potential. Substituting ( |3.7| ) into ( |3.(j ) 

gives the leading-order thermal corrections to the static 



(3.8) 
(3.9) 



and then matching (2^) to 
potential and the singlet normalization: 



where vi°^(r) - V}'''{r) = NcaJ{2r). The power counting of pNRQCD and Eq. (^]^) give 
develops also an imaginary part |^] (see also |^). 



5Zs{r) 



1^ C,a, [Vt'^ {r)) " {v^^ (r) - (r)) , 



vr 
9 



~ afT r 3> mOg. We recall that at higher orders 



3.2 Singlet spin-orbit potential bVi^ssa in pNRQCDjjrpj^ 

In 0], all contributions to SYg^ static and non-static, that contribute to the spectrum up to 
order moc^ were computed. Up to that order, no spin-dependent corrections are relevant. 
The aim of this section is to compute the leading spin-orbit terms in bVs. In particular. 
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Figure 2: Diagrams contributing to SVLS.a- The dot stands for an insertion of the spin-orbit 
potential proportional to Vls sa (left and right diagram) or to Vlsoo (middle diagram), all other 
symbols are as in Fig. |l|. 

we will compute the leading thermal correction, 5Vls sa, to the centre-of-mass momentum 



dependent spin-orbit potential ViSsa, defined in Eq. (2.7). The computation follows the 



same line as the one for we calculate thermal spin-dependent corrections to the 

pNRQCD singlet propagator, and match it to the singlet propagator in pNRQCDjjrpj^. 
We identify the following set of contributions to 5Vls sa ■ 

6VLSsa = 6VLS,a + 6VLS,b + ^3,0 + 6VLS,d + ms,e , (3T0) 

where 

(1) 6VLs,a comes from inserting a spin-orbit potential in the singlet or octet propagators 
of the diagram in Fig. ^ 

(2) 6VLs,b comes from replacing one of the two chromoelectric dipole vertices in Fig. || 
with the chromomagnetic vertex proportional to cpV^^'^^ in Eq. ( p. 13 ) and inserting 
a centre-of-mass kinetic energy into the octet propagator; 

(3) 5Vls,c comes from replacing one of the chromoelectric dipole vertices in Fig. [l| with 
the vertex proportional to CgVs^oa^ in Eq. (|2.15| ); 

(4) dVLs,d comes from replacing one of the chromoelectric dipole vertices in Fig. [l| with 
the vertex proportional to V^^'^} in Eq. (|2.15|) ; 

(5) SVLs,e comes from replacing one of the chromoelectric dipole vertices in Fig. [l| with 
the vertex proportional to cpV^^f^ in Eq. ( |2.13| ) and the other one with the vertex 
given by hi];^^ in Eq. ( |2.14| ). 

By explicit inspection, one sees that diagrams with vertices given by the terms proportional 
to vjo'c \ V^l'!^) and Vl^jl{r) in Eqs. (|1|) and (|2T5D, albeit spin-dependent, do not 
contribute to the spin-orbit potential. 

3.2.1 Evaluation of dViSa 



We evaluate the diagrams in Fig. |2|. As in the previous calculation of the thermal correction 
to the static propagator, we expand the octet propagators for » E — h^^ (see Eq. (|3.5| )). 
The leading contribution comes again from the linear term that we treat by means of the 
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Figure 3: Diagrams contributing to SVLS.b- The square stands for tlie cliromomagnetic vertex 
proportional to CEVso b"^ i''^)^ cross for a centre-of-mass kinetic energy insertion, and all other 
symbols are as in Fig. |l|. 



identity (3^). The left diagram of Fig. ^ gives 



2vr^ 2 i 



E-K 



(0) 



-r{r\E 



hf^}-{v!;'\r)-V^'){r) 



1 (rxP)-(^W-^(2)) ^^ , 1 

X VLSsair)^-^^, 



(3.11) 



the right one gives 



CpasT VLSsa[r) 



9 



E-h 



(0) 



4m2 



i{r^i?-/.W}-(y(°)(r)-yi0)(r) 



1 



(0) ' 



(3.12) 



and the middle one gives 



271" „ 2 



E-h 



(0) 



r(2 



Am? 



VLSoa{r) 



E - hf^ ' 



(3.13) 



where we have kept only terms relevant at order 1/m?. Matching to the pNRQCDjjrpj^ 
propagator ( |3.2D , expanded according to ( |2.6| ), we observe that the terms proportional 
to iVo'^^ — yj^r^ in ( p.llj ) and ( 3.12 ) cancel against one insertion of 5Vs'^\r) and one 
of the spin-orbit potential in the pNRQCDjjrpj^ propagator, while the term proportional 
to {E - /ii°VV2 in ( pTl ) and the one proportional to r'^{E - hf^)/2 in ( |3T2D cancel 
against the term ^5Zs,i/{E - hf^) x [spin-orbit potential] x 1/{E - /ii°^)| in (^^. The 
expression of 5Zs can be read from Eq. ( |3.9D . What is left gives the leading-order thermal 
correction, coming from the diagrams in Fig. |2|, to the centre-of-mass momentum dependent 
spin-orbit potential: 



SVLsAr) = ^C^F«s (vt^\r)fT\^ {VLSoa{r) - VLSsa{r)) 



(3.14) 



According to the power counting of pNRQCD and Eq. (IJ^), we have that ma^ ^ 
^VLS^a (r X P) • (cr(i) - cr(2))/m2 > maf . 
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3.2.2 Evaluation of 5VLs,b 

We evaluate the diagrams in Fig. |3|. Their thermal contribution to reads 



xr 



d'^k i (P-k)2 1 



(2^)4 ^ _ f^m -ko + iv ^"^ E- h^o'^ -ko + iv 

27r(^(A;2-A;2)nB(|A;o|)(a«^-a(2)j), (3.15) 

while the non-thermal part vanishes if regularized in dimensional regularization. The fac- 
tor 2 follows from the fact that the two diagrams give the same contribution at order 
The octet propagators may be expanded according to Eq. ( |3.5D ; considering that, 
besides the two octet propagators, the integral in ( |3.15| ) is odd in ko, the leading non- 
vanishing term coming from their expansion is —2{E — h^^)/{—ko + ir])^. The factor 
E — h^o^ contains a part, E — h^s \ that contributes to the singlet normalization, and a 
part, Vs'^'^ — Vo^\ that contributes to the spin-orbit potential. The octet centre-of-mass 
kinetic energy, (P — k)'^/(4m), contributes to the spin-orbit potential only through the 
term — P-k/(2m). With this in mind, we match ( p.4[ ) to ( |3.2D and obtain the leading-order 
thermal correction, coming from the diagrams in Fig. |3|, to the centre-of-mass momentum 
dependent spin-orbit potential: 

^VLS^r) = -fcFa,CFVj^''\r)v}'J\r)T' (yW(r) - V^ir)) . (3.16) 

Considering that the matching coefficients cp, vjo'^^ and V^^'l^^ are one at leading order, 
the size of the correction is mOg ^ ^VLs,b (r x P) • (cr^^^ — cT'^^~)lrr? ^ maf . 

3.2.3 Evaluation of bVLS.c 

The calculation of 5Vls,c is at this point simple: there are two contributing diagrams, which 
may be constructed by replacing one of the chromoelectric dipole vertices in Fig. |l| with the 
vertex proportional to CsYslia^ in Eq. (|2.15 ). Since this vertex contains a chromoelectric 



field as well, the integration is exactly the same as the one performed in Eq. ( |3.4| ). The 
only change is in the prefactor of the integral. Matching to the pNRQCDjjrp^^ propagator, 
we obtain at leading order 

bVLsAr) = '^^Cpa,c,Vt^\T)V^l^^\T)T^ (v^^\t) - ^^(r)) , (3.17) 
which, considering that c^, Vso'^^ and Vsoa^ are one at leading order, has the same size as 

3.2.4 Evaluation of W^sA 

The calculation of bVi^s d is similar to this last one, but with the vertices proportional to 



CsVsoa replaced by the ones proportional to V^^yi in Eq. ( p. 15 ). The leading-order result 
rGcids 

5V,sAr) = f C^«sKi°'^)(OV;?^^(OT^ {V^'\r) - vP{r)) . (3.18) 
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Figure 4: Diagrams contributing to SVls.c- The diamond stands for the vertex hi],'^\ given in Eq. 



(2.14), and aU other symbols are as in Fig. ^. 



3 



Considering that Vs^o'^^ is one at leading order, but that V^^'^} is at least of order a^, 
6VLs,d{f) is suppressed with respect to 6VLs,a, SyLS,b and 6Vls,c by, at least, a factor a. 

3.2.5 Evaluation of 6VLs,e 

We evaluate the diagrams in Fig. Their thermal contribution to T,s{E) reads 

fPxr)*\ r d^k i 



x(iA:')e-'''^e^"^(-iA:") (^5'' - 2Ti6{kl - k')nB{\kQ\) {a^^^^ - a^^)^) , (3.19) 

while the non-thermal part vanishes if regularized in dimensional regularization. The factor 
2 follows from the fact that the two diagrams give the same contribution at order 1/m^. 
The octet propagators may be expanded according to Eq. (p.5|): the linear term in S — 



h^o^ contains a part, E — h^s \ that contributes to the singlet normalization, and a part, 
Vs^^ — v!f'\ that contributes to the spin-orbit potential. This last contribution reads 

SVLsAr) = Y^Fa,^FV^y\r)V^l';\r)T' (^^(r) - ^^(r)) . (3.20) 

Considering that, according to ( p. 161) , the matching coefficient Fso^'^'* is equal to Vso'^\ 
SyLS.e exactly cancels with SV^s^b in the sum ( 3.10| ). 



3.2.6 Summary 

In summary, the leading thermal correction to the centre-of-mass momentum-dependent 
spin-orbit potential, 

(r X P) • (cr(l) - cr(2)) 

-dependent ^ 2 

-5VLSsa{r), (3.21) 



is the sum of Eqs. (|1|), (pl ), ([sTzl) , (|1|) and (|]20|); it reads: 



SVLSsair) = -^CFa,Vt^\r)T''[vt^\ry{VLSoa{r)-VLSsa{r)) 
[csVt^ir) + V^'Jir)] (F(°)(r) - F (°)(r)) } 



+ 



-CfA^c— + higher orders . (3.22) 
6 r 



In the first equality of ( 3.22| ), the matching coefficients of NRQCD and pNRQCD have 



been kept unexpanded; this amounts at having provided an expression for the spin-orbit 
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potential that resums contributions from the scales m and mas, while it is of leading order 
in the temperature. In the second equality, we have kept only the leading terms in the 
NRQCD and pNRQCD matching coefficients. We note that the contribution coming from 
the term proportional to V^'^'i^ii is negligible, of the same size or smaller than subleading 
thermal corrections that we have neglected throughout the paper. 



4. Gromes relation at finite temperature 

After having computed the leading contributions to 5Vls sa, we can now check whether 



these new terms fulfill the Gromes relation ( 2.10 ) or not. This corresponds to verifying the 
equality 



mssa{r) = (4.1) 



We use the expression of SViSsa provided by the first equality in Eq. ( |3.22|) that keeps 
unexpanded the matching coefficients of NRQCD and pNRQCD. If we make use of the 
relations ( |2.10| ) and ( p.l9 ), which are exact, then 5Vls sa may be rewritten as 



OVLSsa[r) = 

2r 

+ |c^a.yi°'^)(r)T2 (c.yi„Y)(r) + yi°'«(r)) (^^(r) - V^'){r)) ,(4.2) 
which shows that the Gromes relation is violated by an amount, which at leading order is 
9 r 

4.1 The spin-orbit potential in a covariant model 

In order to understand the origin of the observed violation of Poincare invariance, it is 
useful to consider at zero temperature the case of a massive gluon, whose mass, rUg, is 
such that mos ^ nig » mOg. The massive gluon contributes to the potential, but clearly 
it does not break Poincare invariance. To see this let us evaluate the corrections to the 
spin-orbit potential. The diagrams contributing to the spin-orbit potential are the same 
of those considered in the previous section, only now the gluon propagator reads (in the 
unitary gauge) 

V - ^-^) ■ (4.3) 



k^ — kp' — nig + IT] \^^^ y 
The contributions to the static potential, 5Vs^\ and SVis^a, ^Vls^c and 5VLs,d depend 



on the correlator of two chromoelectric fields. They are proportional to 

h)o' - ko + irj '■"(J 



/ WfE-h^'^'-k. + ^r^kl-k^-ml + ^r, " '° " ' ^^'^^ 



where we have regularized the integral in dimensional regularization. Integrating over the 
momentum region k^^k ^ nig means that we are expanding the octet propagator according 
to Eq. ( |3.5D . We will focus here only on the term linear vn E — h^^ , since this is the relevant 
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term in the finite temperature case analyzed in the rest of the paper. ^ It turns out that 
the hnear term vanishes in dimensional regularization (see The reason is that the 

contribution coming from the spatial components of the gluon propagator (proportional 
to {d — 1) /cq ill the equation) cancels against the contribution coming from the temporal 
component (proportional to fc^ in the equation). This is in sharp contrast with the finite 
temperature case, where the term linear in E — h^^ does not vanish (in Coulomb gauge, 
this is due to the fact that only the spatial components of the gluon propagator get thermal 
contributions) and eventually generate a finite thermal contribution to SVj'^\ SyLS,a, SyLS,c 
and dVLs,d (see Eqs. (U), (|3l4|) , (IsTtI ) and (|1|) respectively). 



The contributions to 6VLs,b and 6VLs,e depend on the spatial components of the gluon 
propagator only. Both (5Vl5,6 and SVLs,e get finite contributions from the massive gluon 
but the sum of the two terms linear in E — h^^ vanishes: the same happens in the finite 
temperature case discussed in the previous section. 

The massive gluon example provides a simple case where Poincare invariance is not 
broken. The Gromes relation is trivially realized for terms that are linear in the energy: 
such terms vanish for both the static and the spin-orbit potentials. In the finite temperature 
case, diagrams that depend on the correlator of two chromoelectric fields, like the one shown 
in Fig. m, do not vanish. This is a direct consequence of the fact that the thermal bath 
affects in a non-covariant way the gluon propagator. 



5. Singlet spin-orbit potential SVLSsb in pNRQCDjjrpL 

In this section, we calculate the leading thermal corrections to the spin-orbit potential 
SVls sb: which is the spin-orbit potential experienced by the quarkonium when at rest with 
respect to the laboratory reference frame (we recall that, in our setup, this is also the 
reference frame of the thermal bath). This potential, even at zero temperature, is not 
constrained by Poincare invariance. 

In order to calculate SVissb, we need to consider two new terms contributing to /i^o 
in the pNRQCD Lagrangian: the term 

_f^^(i).^^(9^gB), (5.1) 
4m 

and the term 



^^,cT^'^-[px,gE], (5.2) 

where, for simplicity, we have put to their tree-level values the pNRQCD matching coeffi- 
cients. 

There are three classes of diagrams that contribute: 

SVls sb = SVls (i) + SVls (ii) + S^LS (m) • (5-3) 



^We have explicitly checked that also the leading contribution to the spin-orbit potential, proportional 
to rrig/m^, satisfies the Gromes relation [psl. 
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(1) The first class is similar to the one shown in Fig. |2|, but now the dots stand for 
insertions of the spin-orbit potential proportional to VLSsb (left and right diagram) 
or to VLSob (middle diagram). Vissb and VlSo6 have been defined in Eqs. ( ^ ) and 
(|2.8D , and given at leading order in Eq. (|2.9D . The result reads 



2tt 

SyLS(i){r) = —CFa,T\^ {VLSobir) - V^s sb{r)) . (5.4) 



(2) The second class is similar to the one shown in Fig. ^ but now the squares stand for 
the vertex induced by ( |5.1| ) and the cross for a kinetic energy insertion, p^/m. The 
result reads 

5VLSi^i){r) = -^Cpa^cpT^ (yW(r) - vP{r)) . (5.5) 



(3) Finally, the third class of diagrams is similar to the ones evaluated in Sec. 3.2.3, but 
with the vertex proportional to CgVsoa^ in Eq. ( |2.15 ) replaced by the vertex induced 
by (15.21). The result reads 

5VLSi^^^){r) = yC^as^T^ (yW(r) - vP{r)) . (5.6) 
Summing up all three contributions we obtain 
SVLSsb{r) = yCi^asT^ [r2 {VLSob{r) - VLSsb{r)) - (Fi°)(r) - ^/"^(r)) 

= CpNc—T^ + higher orders , (5.7) 

18 r 

where, in the first equality, we have used that 2cf — — 1 = 0. 
6. Conclusions 

We have calculated the leading-order thermal corrections to the quarkonium spin-orbit 
potentials. These corrections go quadratically with the temperature and are proportional 
to alT'^/r. 

At zero temperature, the spin-orbit potential that depends on the centre-of-mass mo- 
mentum is protected by Poincare invariance. We have computed its leading thermal cor- 
rection in Eq. ( 3.22| ). In Eq. ([4.2| ), this correction has been shown to violate Poincare 



invariance. This implies that order a^T^/r corrections to the quarkonium potential will be 
experienced by the system differently in different reference frames, and, in particular, in a 
frame where the thermal bath is not at rest. 

We have also computed the leading thermal correction to the spin-orbit potential of a 
quarkonium at rest with respect to the laboratory reference frame. Its expression is in Eq. 
( |5.7| ). The potential contributes to the spin-orbit splittings of the quarkonium levels. The 
thermal correction being negative implies a weakening of the spin-orbit interaction in the 
medium. 

In the special case of a bound state made of two heavy quarks with different masses, 
mi S> m2 (like the Be system), the complete spin-dependent potential, up to corrections 
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of relative order m2/mi, is given by {VLSsa + ^L5s6)(r x P2) • ''■2/(4rn.2). Therefore, the 
leading contribution to the hyperfine splittings reads 

36 m2 

where n, j and I are the principal, orbital and total angular momentum quantum numbers. 
The thermal correction has opposite sign with respect to the splitting at zero temperature, 
which implies smaller hyperfine splittings in the medium. 
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